ON GENERALIZED HILBERT MATRICES 



RUIMING ZHANG 

Abstract. In this note, we present a systematic method to explicitly compute 
the determinants and inverses for some generalized Hilbert matrices associated 
with orthogonal systems with explicit representations. We expressed the de- 
terminant, the inverse and a lower bound for the smallest eigenvalue of such 
matrix in terms of the orthogonal system. 



1. Introduction 

The Hilbert matrices are the moment matrices associated with Legendre poly- 
nomials. The generalized Hilbert matrices are the generalized moment matrices 
associated with certain orthogonal systems. It may be interesting to find the exact 
formulas for their determinants and inverses. In this note we provide a systematic 
method from the theory of orthogonal polynomials to find these formulas through 
the explicit representations of their related orthogonal systems. We will demon- 
strate that, once the orthogonal system is explicitly known, then the inverse and 
determinant of each generalized moment matrix are known explicitly. Furthermore, 
we also know a lower bound for the smallest eigenvalue. Since an orthogonal system 
may be derived in many different ways, this method could be very handy. We also 
present four examples to show how to apply this method in various situations. 



2. Main Results 

Let E he a complex inner product space with a sequence of linearly independent 
vectors {wnj^p- ■^'-'^ each nonnegative integer n, the following matrix is positive 
definite, 

where 

There is a unique orthonormal system {pfcl^^o '^ith p„ having positive leading co- 
efficient in Un, which could be found through the Gram-Schmidt orthogonalization 
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process. Each Pn is given explicitly by [2] 
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i/det Gn det G„- 



:det 
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— l,n 



J=0 



det G„ = ]J ( 



det Gn-i 
det Gn 



Let {wnjj^o another sequence of vectors in E related to {u„}^q in the followin 
way, 

n 

then, they are also linearly independent. Let us define 



then 



where G* is the Hermitian conjugate of G„ and 

G„ = i'^j,k)j f.^Q , 



with 
Thus, 

Let 
and 



Cj,fe 



0, k> j. 

n 

det Hn — det G„ Cjj. 



Pn = ^ bn,kVk, 



k=0 



where we follow the same convention as above, 

aj,k = foj.fe ==0, /c > j. 
Evidently, both An and i?n are invertible. Let 
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then, 



and 



Thus, 



e=o m=0 



m=0 



which gives 



We summarize our above discussion as the following theorem: 

Theorem 1. For each nonnegative integer Tij dsstime that Gn; -^nj and 
Cn as defined above, then, 

Hn = GnC*, ^ = B*An. 



Ci. 



3,3 
1 



det G„ = J]^ a^. J, det C'n = J]^ Cj,^-, det -f^n = 

where j4* is the Hermitian conjugate of An . 
Notice that 

■■ = (7j,fe)",fe=o 

with 

n 

l3,k = ^ bi,jae,k- 

£—niaK(j.k) 

From the expressions of det Gn and jj^k , we observe that in the case Uk = Vk, the 
requirement that ak,k > could be disregarded in the actual computations, for, if 
we replace the pair {uk,Vk} by {€kUk,€kVk} with |efe| = 1 we won't change Hn- In 
some situations, the entries for G„ are too complicated to be useful. This is the 
reason we won't compute G„ for the Askey- Wilson polynomials, we compute Hn 
instead. But G„ could be recovered easily via 

Gn = Hn (C*) , detG„ = Yfn — — ■ 

lij=0 ^3,3 

Let As be the smallest eigenvalue of G„ , then 4- is the largest eigenvalue of 
Gn^ . Since G~^ is positive definite, we have 



^ -Wg-H < 



J2 i^'-j'i^ 



\ h3=0 



with 
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and 



Hence, 



and 



^ n n 



Another lower bound could be found by considering the 1 1 • 1 1 cxd ■ Let {xo ,.. .,Xn) 
be an eigenvector of corresponding to A^^, then we have 



f ^0 \ 




/ Po,o ■ 


PO,n 






\ Pn,0 ■ 


Pn,n 







1 


\ Xn J 



/A \ 

> ,PO,kXk 
fe=0 



^ ^j Pn,kXk 
\ fe=0 J 



then, 



/ Xo 



( ^ \ 

2_^P0,kXk 
k=0 



^ ^j Pn,kXk 
\ fe=0 / 



0<j<n 











U=o ) 


\ Xn J 



Observe that 



n n 



^\Pj,k\ < ^^\ae,j\ \ai,k\ = ^\ae,j\^\ai,k\ > 



fe=0 



fe=0 ^=0 



^=0 fe=0 
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and 



0<j<n 



. e=0 k=0 



< 



{i«^ji}y]i«^,fci 

0<j<n '—^ 



£=0 



fc=0 



n I n 



< X! XI I'^'^ji 'X! I'^'^^'^i 



e=o I j=o 



k=0 
2 



n ( i 

Theorem 2. For each nonnegative integer n, assume that Gn and An, as defined 
above, then. 



Xs > max < 



1 



1 



E"=oE,toKir'ELo{E-=ok^,.l}\ 

where As is the smallest eigenvalue of Gn ■ 

The second lower bound is particular interesting when the generalized orthogonal 
system are certain orthogonal polynomials, since it could be expressed in terms of 
the orthonormal polynomials. 



Corollary 3. Let 



Pn{x) 



^an,kx'', n = 0, 1,... 

k=0 



be the orthonormal polynomials and rUn be the ii-th power moment with respect to 
a probability measure d^{x). If there is a complex number zq with \zo\ = 1 such 
that all of 

<^n,fc'2o ) fc = 0, 1, . . . 
are of the same sign, then the smallest eigenvalue Ag of the matrix 



has a lower bond 



Gn — {^i+j, 



i,J=0 



ELobK^o)p- 



The first observation for the last corollary is that it seems such zo always exits. 
For the symmetric orthogonal polynomials, from the three term recurrence we see 
that zo = i. But we don't know how to prove this for general cases. The second 



ON GENERALIZED HILBERT MATRICES 



6 



observation is that if all the polynomials are real, then the special case of Christoffel- 
Darboux formula gives 



m=0 

which gives us 

A, > 



{Pn+l{zo)Pn {Zo) - K(zo)Pn+l (^o) } 



ari+l,Ti+l 



{Pn+li^o)Pnizo) - p'nizo)Pn+lizo)} ' 

and this may be useful to find the asymptotic behaviour of the lower bound. 

2.1. Matrices associated Muntz systems. 
Theorem 4. For ti G N and {oq, ai, . . . } C C, the matrix 

1 



(2.1) 

has determinant 



det 



aj +ak + 1 



j,k=a 



a, + Qffc + 1 



n n 



n 



fc-i 

Y[\ak - aj\ 



j.k=0 k=Q 



(l + 23fi(afe))[]|afe + cE-+l| 



3=0 



Under the condition 

aj ^ Q!fe, 3?(aj) + ^{ak) ^ A; = 0, 1, . . .n, 

the matrix (|2.ip is invertible, and its inverse matrix {jj,k)"k=o element 



m— max(j,/c) 



(1 + 25R(a„)) Yl (aJ + ar + l){ak + a; + 1) 

r=0 

m — 1 m — 1 

n - ^) n ("'^ " 

p=0 Q = 

p^j q^k 
When the matrix l|2.ip is positive definite, its smallest eigenvalue has a lower hound 



K > < 



E 

i=0 



3=0 



fc = 



More generaHzed matrices associated with a generalized Miintz system: 
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Theorem 5. Given a,b,c ^ R and distinct complex numbers {a/c}^o ' ^^^^ 
positive integer n, the matrix 



(2.2) 
has 



1 



coijOik - a{aj + Uk) - bj .^^^ 



det 



XjUk - a{aj + ttfe) - b J J j^^g 



fc-i 



(a^ + 6c)^^^ Jl (c|afc|2 - 2a3?(afe) - b) Y[\ak - 



3=0 



Under the conditions 



+ bc^ 0, cajak — ci{aj + ak) — b ^ 0, aj ^ ak 

for j,k = 0, ...,rt, i/ie matrix l|2.2p is invertible, and its inverse {'Jj,k)"k=() 
element 

m—l 

^ (cjamp — 2adi{am) — b) {carCtJ — a{ar + aj) — b){ca^ak — a{a^ + ak) — b) 

r=0 



p=0 q^O 
pj^j p^k 

When the matrix (|2.2p is positive definite, its smallest eigenvalue has a lower 
bound 

2x -1 



A, > < 



E 

1=0 



^ ^cja^p - 2a3fJ(a£) - 6j_ J_|cQ;jQ;fc - a{aj + ak) - b\ 
i^o 

£-1 



3=0 



{a^ + bcy/^ n 

= 

2.2. Matrices Associated with g-Orthogonal polynomials. Recall that for 

a e C and g e (0, 1), [2J, 



/ N (a; 9)00 
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j 



-, 0<j<m, 



and 



(ai,a2, ...,an;q)m = JJ(afc;g)m, m e Z,n e N 



k=l 



for ai, 02, a„ e C. 
For J/ > -1, let 



< j,y,i{q) < ivaiq) < ■ ■ ■ < iuAq) < 

-V t(2) 



be the positive zeros of z '^Jl^'{z;q), where the Jackson's g-Bessel function J^\z;q) 
is defined as 



Let us define 



where 



wi 



Sn u —7 .2„ / \ ) 



Wl 



M+1 (>,^(g);g) 

d.jl^^ {z;q) 



\z=ju,iiq) ' 



we have the following result: 

Theorem 6. For v > —\ and n G N, the matrix 

(2-3) {sj+k+i,uTj^k=a 
has determinant 

2-n(n+l)^n(n+l)(4ri+6iy+5) 



det {Sj+k+l,u)j^k=0 



rn—1 

and its inverse matrix (7j,fc)J;.^o element 



£— max(j,fc) 

^(2j+2/c+l)£ (g-+l;g),_J 







' e + k' 




9 


i-k 



The smallest eigenvalue of the matrix (|2.3p /las a lower bound 



«=0 
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where 



^ (-l)"(9^g;9)«-fe(g;'7)2fc 



Matrices associated with Askey- Wilson orthogonal polynomials: 

Theorem 7. For n G N, i/ie matrix 

(a;g)j+fc V' 
(a/3; g)j+fc/j-fc^o 
/las determinant 



(2.4) 



det 



(a/3; (j)j+fc 



„+l)/2^„(„^-l)/3-Q(g^^^^. 



]J(a/3g™-i;g)„.(a/3;g) 



2m 



m— 1 



Under the conditions 



a ^ 0, q ^ 0, q, a, P ^ q , k — 0, . . .n, 
the matrix (|2.4p is invertible, and its inverse matrix ijj,k)"k=o element 



m 




m 


3 


q 


k 



(agJ+fc)™(g, /3; q)„^{af3; g)2m(a/3(?™-i; g). 



_ ^^ly+^qm + m " - {aPq"^-';qUa(3q"^-\qh 

^'"^ (a;g)j(a;g)fe 

When < g, a, /3 < 1 i/ie matrix \2A\ is positive definite, its smallest eigenvalue 
of the matrix has a lower hound 

-1 



where 



K > 



(«/3;g)c 
(a;'?)oc 



E 

.£=0 



p^,i-l;aq-\Pq-'\q) 
hi{aq-^ , I3q-^\q) 



Pn[x;aq , (3q |g) = 20i ( ^ ; g; 



/i„(ag \/3g ^\q) 



(a/3;g)oo (l-a/3g ^)(g, /3; g)„a" 
(a;g)oo (1 - a/3g2»-i)(a, a/3g-i; q). 



Remark 8. One could get their classical counterparts for the formulas in Theorem 
[6] and Theorem [7| by passing the limit g — > 1~ with proper normalizations. 



3. Proofs 

In this section we prove our results under the same restrictions for the orthog- 
onal polynomials. But the results clearly hold for more general cases since all 
the expressions involved except the third example are rational functions in their 
parameters. 
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3.1. Proof for Theorem[4l Given a sequence of distinct complex numbers {ctk}^^Q 
satisfying 5R(afe) > — ^, the orthogonal Miintz-Legendre polynomials are defined as 

m 



n-l 



Ln{x;ao, . . . , q;„) = ^ Cn,kX°"' , 



/c=0 



n ("'^ ~ 



for n e N and 

Lo(a;;ao, . . . ,a„) = 

They satisfy the following orthogonal relation 
1 



Ln{x] ao, ■ • ■ , an)Lm{x\ ao, . . . , am)dx = 



/o l + 2>H(a„j 

for n, m e {0} U N. We take 

and 

The orthonormal polynomials are given by 



Pn{x) =v^l + 2^{an)Ln{x] ttQ, . . . , a„). 

Then, 

Tl-l 

On,*; — OnM — — [ 

n 

J = 

Theorem |4] follows from Theorem [T] and Theorem [21 

3.2. Proof for Theorem [5l Assume that 

+ hc> 0, 

and 

CQfjOfc — a{aj + Sfe) — 6 > 0, 
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the following orthogonal Miintz polynomials exist, jSj 

n-l 

n( 



n-l 

ttfc =i 



J=0 



n ~ 

They satisfy the orthogonal relation 
where 

(a2 + 6c)" 

(c|a„|2 - 2a3fiK) - 6) [] |mj - a|2 

i=o 

Take 

From [3] we have 

^''^ cajOfe - a(Q;j + Ok) - 6 
The orthonormal polynomials are 

n 

Pn{x) = y^fln.fca:"'', 

where 

0^71 k Oil i; 



Then, Theorem [5] follows from Theorem [1] and Theorem [H 

3.3. Proof for Theorem [6l The even g-Lommel polynomials satisfy the orthog- 
onal relation [2] 



Let 



- — rT;? - 



„2ni/+n(2n+l) r 



then 

Ij.k = Sj+k+l,vi Cj,k — Sj^k- 

The orthonormal polynomials are given by 



^1 _ g2n+i.+l 
PA^) = ^n.+n(n+l/2) fe2r..+l (a:; g) . 
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hence 

Then, Theorem [6] follows from Theorem [1] and Theorem [2l 

3.4. Proof for Theorem O We present two proofs to this example. Our first 
proof uses the Askey- Wilson polynomials, while the second uses the little q-Jacobi 
polynomials. For some mysterious reasons, these polynomials yield essentially the 
same matrices. 

For each nonnegative integer n, the Askey- Wilson polynomial has the following 
series representation 

a„(x; t\q) = ti'^{tit2, hts, tiU; q)n 

X4'?^3 i + i + + i \q,q), cose, 

where the basic hypergeometric function r</>s with complex parameters ai, ...,ar',bi, ...^bs 
is formally defined as, [2J 



^ {q,bi,...,bs]q)n V / 



It is well known that an{x]t\q) is symmetric in the real parameters ti,t2,t3,t4. 
Under the condition max { | ti |, |t2|, 1^3 1, 1^41} < 1, the Askey- Wilson polynomials 
satisfy the following orthogonal relation 

nl 

a^{x;t\q)an{x;t\q)w{x;t\q)dx = /i„5^,„. 



where 



and 



('z"+^'?)oo n (*j*fe'?";9)oo 

l<j<k<4 



w{x;t\q) ^ — --j==, x = cost 



Let 



it is known that f2] 

(be , be , gj„ _ {q , ge , ae j^q 
{ab,b/a;q)n ^ {q,ab,q^-"'a/b;q)k 



for a ■ b ^ 0, hence, 

k 

n /. , h \ I T^2 \ I ^2 ^ 



n 

k 



1 



and 
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or 



li.k = 



Ij,k = j w{x;txq\t2q'^ ,tz,ti;q)dx, 



^'^ [q, tit2q^+'^, tit^qi, tiUqi ,t2hq'^ ,t2Uq'^ ^hU; q)oo ' 
The orthonormal polynomials are 



aj{x;t\q) 



with 



and 



Then 



Pj{x) = 



_ {tit2,tit3,tit4;q)„iq ",^1^2*3*49" ^■,q)kq'' 



(ti)" V7i^(g, tit2,tit3, tiU; q)k 

{t2ti,t2h, t2t4; q)n{q^'\ tit2tzUq'^~^]q)kq^ 



{t2)"'\/K{q, t2tl, t2H, t2U\ q)k 



(3.1) det 



A.Tj 27r (titag'"-^)" {tit2hUq^"' ; q)^ 
aetii„= — , 

l<j<k<4 

which could be simplified to 

{tit2;q)j+k Y 
Stit2t3t4;q)j+kJj^k=o 

_ A (tlt2g'"~^)'" (g, tit2,t3t4; q)m 
m=0 (*l*2*3t4; q)2m{tlt2t3Uq"'~'^; q)m ' 

For any n = 0, 1, let e, cq, Ci, c„, do, di, rf„ be non-zero numbers and, 
^ = i^jk)lk=0 ' ^ = iyjk)lk=0 ' ^ = ihk)lk=0 ■ 

If 

XY = 1, 

then, 

XY = 1, 

where 

X = {eXjkCjdk}2k^Q , 

and 



edjCk J j 
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Using the above trick, we simplify the inverse pairs down to 

fti^2; q)j+k 



and its inverse matrix {i^j.kYjk^o ^^^^ element 



(3.2) 



{tlt2;q)j{tlt2]q)k fr'n ^3^4; g)r 



m 




m 


j 


1 


k 



m=0 

{tit2t3Uq"'-^;q),{ht2t3Uq"'-^;q)k 
itit2q^+''r {tlt2UU- q)2m{tlt2UUq'''~'^-q)r 

Except the lower bound, the assertions of Theorem [7] follows from l|3.ip and p. 2 
by the change of variables, 

a = tit2, P = ^3^4- 

Assume that 

p-i{x; a, b\q) = 0, po{x; a, b\q) = 1, 
the little (?-Jacobi polynomials {Pn{x; a, b\q)}^^Q have the orthogonal relation [2] 



E ^^^7^M^^'™('?'; b\q)pn{q'';a, b\q) = /i„(a, b\q)5„ 



k=0 



{q;q)k 



for m,n> The moments are given by the formula 

{bq;qU{aqrq" 



or 



m=0 



^J■n 



iq;q)rr. 

{abq^+^-q)r_ 



(aq^+^;q)^ ' 

by using the g-binomial theorem |2|. The orthonormal polynomial 

{-l)"Pn{x; aq, bq\q) 



Hence, 



Pn{x) 



\/Ma7%) 



^Jhn{a,b\q)[q,aq] q)k 
Theorem [7] follows from Theorem [T] and Theorem [2] with a change of variables 

a — aq, /3 = bq. 
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